Quantum description of light pulse scattering on a single atom in waveguides 
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We present a time dependent quantum calculation of the scattering of a few-photon pulse on a 
single atom. The photon wave packet is assumed to propagate in a transversely strongly confined 
geometry, which ensures strong atom-light coupling and allows a quasi ID treatment. The amplitude 
and phase of the transmitted, reflected and transversely scattered part of the wave packet strongly 
depend on the pulse length (bandwidth) and energy. For a transverse mode size of the order of 
A 2 , we find nonlinear behavior for a few photons already, or even for a single photon. In a second 
step we study the collision of two such wave packets at the atomic site and find striking differences 
£S] between Fock state and coherent state wave packets of the same photon number. 
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I. INTRODUCTION 

Identifying and realizing systems with strong coupling between light and matter is one of the central objectives of 
' current research in quantum optics. In the strong coupling regime the coherent interaction between a few atoms and 
, the radiation held of few photons takes place on a fast time scale and is not masked by the dissipative coupling to the 
environment. A weak quantum field can change the internal quantum state of an atom and, vice versa, a single atom 
is able to have an appreciable influence on the light field. Both effects have, evidently, many interesting practical 
i^: 1 implications including single atom detectors as well as nonlinear micro-optical elements on the single-photon level, 
\ which are in the heart of all-optical quantum information processing schemes jjj, e.g. to facilitate Bell state detectors 
■ for photons 

So far optical experiments in the strong coupling regime have been performed with atoms put in tiny high-Q cavities 
H 0]. Here the origin of the strong coupling is the small volume of the discrete cavity modes. The local electric 
field sustained by a single photon becomes then very large. Strong coupling in the microwave regime yielded several 
spectacular results q{. However, the electric field enclosed in three dimensions allows only limited accessibility, 
which complicates direct detection of the light and limits the time scales. In addition, the scalability to a large number 
^ \ of such elements seems rather difficult. 

Hence, there is still a high demand for alternative setups which allow controlled manipulation and strong (nonlinear) 
coupling of single photons. It seems natural to try to lift the confinement of the light field in one particular dimension 
and consider the interaction of transversely confined propagating fields with nonlinear optical elements. If a photon 
forms a very short wave packet along its propagating direction, we again get a small total volume and a high field 
per photon. Conceptually the most simple nonlinear element for such wave packets is of course a single atom with a 
resonant dipole transition. An atomic dipole can be associated with an effective radiative cross section of a a — 3\ 2 /2ir. 
Hence, as a first guess one would expect that it is sufficient to simply focus a pulsed light beam down to a cross section 
of the same order of magnitude as a a to enter the strong coupling regime. This is not too restrictive as beams can 
be focused down to a spot size smaller than A 2 . Unfortunately, this line of argumentation is too naive and a precise 
calculation for such a strongly focused beam shows that the light will propagate through the atom without being 
noticeably influenced [Q . The reason for this can be traced back to the large range of transverse wave vectors involved 
in the dynamics yielding rapid variation of input phase and polarization. Another way to enhanced effective atomic 
cross section could, in principle, consist in using a quantum degenerate gas as an optically dense medium ||. Here 
the use of many atoms leads to quite stringent practical limitations. 

It is the advent of highly developed nanofabrication techniques that may open the way to new geometries avoiding 
these problems and facilitating setups to implement quasi one-dimensional scattering processes on dipoles. Minia- 
turized waveguides on surfaces or optical fibers can be fabricated with cross sections on the order of the optical 
wavelength square. Photons traveling within are well confined in the transverse direction, while the field remains 
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uniform in the longitudinal direction of propagation. The problems of the beam divergence and rapidly varying wave 
vector, encountered in free space, are missing in the waveguide setup. Hence one can hope that a single atom placed 
in the field of such a mode will have considerable influence on the dynamics. The investigation of the prospects of 
such a setup is the subject of the present paper. 

We will study the time dependent interaction of an atomic point dipole modeled by a two-level atom and the 
quantized electromagnetic field in a single transverse mode waveguide. We use the term "waveguide" in a general 
sense, avoiding any precise geometric specification. Instead, we set up a model that accounts for the generic features of 
these devices. The central property is that the atom interacts with a continuum of transversally confined propagating 
modes. A similar system has been studied to investigate the influence of input photon-statistics on single-atom 
absorption || . For this the theory of cascaded open systems JlC , 1 1 as an extended form of the input-output formalism 
was applied. In contrast to the stationary scattering scenario, we will concentrate here on the time-dependence of 
the scattering process, i.e. we consider wave packets consisting of one or a few photons impinging on an atom. This 
allows to find the dependence of the atomic dynamics and the amount of scattering on the light pulse bandwidth and 
energy beyond the narrow-bandwidth limit considered in Rcf. ^ . For very short pulses the energy of a single photon 
is strongly concentrated in space and thus saturation is expected to prominently influence the scattering process even 
for weak energy incident wave packets. As a consequence, we have to take the full quantum dynamics of the dipole 
into account, including saturation to all orders. 

To treat this problem of wave packet scattering we develop a method based on a set of Heisenberg-Langevin 
equations for the field and atomic operators. These equations, relying on first principles of quantum mechanics, 
will be transformed into a simplified form by means of a Markoff approximation. This approximation closely relates 
our approach to the cascaded open systems method. Nevertheless we are primarily interested in the field dynamics, 
called "channel" in the scheme of cascaded systems. Hence we will keep the waveguide electric field as an explicit 
dynamic component of our model. Let us mention here that some numerical examples calculating the electric field of 
the propagation of a single-photon wave packet through a single atom have been obtained before Similarly the 
dynamics of a 2D 1-photon wave packet and many two-state atoms modeling a beam splitter was numerically solved 
& 

Nevertheless, the systematic and exhaustive discussion of the parameter space by a purely numerical approach seems 
hopeless. Fortunately in our method, we are able to get explicit analytical expressions for the scattered waveguide 
field. This result exhibits the scaling with the ratio of the transverse mode size and the atomic cross section. Thus we 
can easily reveal how this ratio governs the energy redistribution in a scattering process. In addition, physical insight 
in the phase properties and the pulse shape deformation can also be gained. 

As a possible application of our model we address the problem of atom-mediated nonlinear photon-photon inter- 
action. So far, significant coupling has been reached with the help of a high-Q cavity |Q or might be reached by 
increasing the number of atoms jg, [l5| . The waveguide-atom interaction could be the basis of nonlinear pulse am- 
plifiers or even that of photonic quantum gates. Here we limit the study to a prototype interference experiment in 
which two wave packets impinge simultaneously on the same atom. 

The paper is organized as follows. In Section II we present the model and introduce the basic parameter describing 
the coupling of the atomic dipole and the waveguide, and finally, calculate the scattered electric field generally within 
a Markoff approximation. In Section III the transmitted and reflected Poynting vectors are studied for various initial 
states of the wave packet. These states include coherent states and the single-photon Fock state. We study the role 
of the pulse bandwidth in the scattering process and the saturation nonlinearity for multi-photon wave packets. The 
phase properties and the pulse shape deformation is also discussed in this section. Then, Section IV is devoted to the 
interaction of light pulses. We conclude finally in Section V. 



II. THE MODEL 



Let us consider a single branch of modes propagating in the +z direction of a lossless waveguide. The longitudinal 
wave number k of the modes is assumed to obey the simple dispersion relation k — co/c, which is valid far from the 
branch threshold |L6| . The effective transverse cross section of the modes 

= J dxdy\f k (x,y)\ 2 

\fk{x ai Va)\ 2 

where fk(x, y) is the transverse mode function and (x ai y a ) the position of the atom, is approximately constant in the 
range of the relevant longitudinal wave numbers. The details of the transverse mode profile do not play an important 
role in the forthcoming calculations as long as the atom may be treated as point-like. Only the value of the mode 
function at the position of the atom enters in the calculation. This value has been incorporated in the definition 
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of the effective cross section A. Similarly, the polarization properties of the field can also be omitted provided the 
field is closely uniform across the spatial wave function of the atom. No assumption is made on the atomic position 
relative to the waveguide. In principle the atom can sit inside a hollow waveguide, it can be embedded in the dielectric 
material of a fiber or can be even outside a dielectric interacting with the evanescent field close to the surface. 

The continuum field quantization follows the theory presented in Ref. |17| . The electric radiation field is decomposed 
into positive and negative frequency parts 

E(z,t) = E ( -+\z,t) (z,t) , (1) 

where 



(z, t)=ij"° ^4^4 («"W e ~ Mt ~ Z/C) + Ut)e- Mt+Z/C) ) 



(2) 



We separated the modes of the two counter-propagating directions into sets of and b u , the first corresponding to 
the forward, the latter to the backward propagating modes. The electric field is given in the interaction picture. The 
field amplitude variables a u (t) and (t) describe then the time variation due to the interaction with the atom. They 
follow the usual bosonic commutation rules, 

a4t),al(t)] = [b u (t),bl(t)] =%-«/) , (3) 

all the other commutators vanish. 

The atomic dipole, again in interaction picture, reads 



leg 



(a~e-' luJAt +a+e luJAt ) , (4) 



where the operators a± together with <j z = (<r + a- — er_<7 + )/2 form a pseudospin obeying the spin-i algebra. We 
assume that the dipole moment is oriented parallel to the field polarization at the atomic position yielding maximum 
coupling. 

The dipole interaction Hamiltonian in the rotating wave approximation is given by 

Hj = -ih J dujg u (a+{a w j" ZA / c + 6 we -*"Wc) e -<(«--«-A)t _ h ^ ; (5) 
where za is the position of the atom, and the coupling constants are 



9 - = v ^kx de9 ■ (6) 

Note that the dimension of the coupling constant g^ is not a frequency but 1 /y'sec. The use of the rotating wave 
approximation is justified as we consider light pulses with bandwidth much smaller than the frequency uja- 
The evolution of the system variables is governed by a set of coupled Heisenberg-Langevin equations 

= s^-e-^^e^-^) 1 , (7a) 

at 

= s^.e^/V^-^* , (7b) 

at 
d 

-r (J - 

dt 

d 

dt z 

Besides the terms originating from the interaction Hamiltonian (^|), we account for the interaction of the atom with an 
environment. This coupling results in a dissipation process with decay rate 70 and with associated noise represented by 
the operators £. The physical role of this 70 decay process is that it provides a channel for the transverse scattering, 
i.e., when photons are scattered by the atom out of the waveguide. The environment consists of the free space 
radiation modes H]. Of course, the presence of the ID waveguide slightly perturbs the surrounding mode structure 
with respect to the simple free space case. This effect, which depends on the specific choice of the waveguide geometry, 
is neglected and, for the sake of simplicity, in the numerical examples we will use the free-space spontaneous decay 
rate 70 = uJ A dl g /(6neohc 3 ). 



- 7o a_ + 2ct z J dLog^a^ ZA ' c + b ^ e -^ A /c^ e -iW-w A )t + |_ ( (7c) 
-2 70 ( ( 7 Z + 1/2) - / dwg u (a+(a^ ZA / c + ^e^i/yMiX + h .c) + |, . (7d) 
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The waveguide modes form a one-dimensional continuum. The above equations can be transformed then into a 
much simpler form by means of a Markoff approximation. On integrating Eq. (Qa), the waveguide field amplitudes 
arise in a form decomposed into a free field part and into a part radiated by the atom 111] , 

M*) = aM+9^-^ ZA / c f dt'a-it'Y^-^' , (8) 



J to 

respectively. The back propagating modes are decomposed analogously. The Markoff approximation is invoked to 
describe the back action of the second term on the atom. It can be applied because the continuum is broadband 
around the atomic frequency. The free field term is usually identified with a Langevin-type noise source. In the 
waveguide we can not make this step since we will consider initial field states different from the vacuum. Altogether, 
the atom is subject to the effect of the free field and to a relaxation into the waveguide continuum. The polarization 
operator cr_, for example, follows 

= -(70 + fi)<T- + 2ct z J du;gMto)e WZA/c + b^e-^'^e^"-"^* + £_ , (9) 

where a new damping rate 71 appears, and the second term contains the free field contributions. The vacuum 
frequency shift induced by the waveguide modes, accompanying the relaxation, is assumed to be already incorporated 
in a rcnormalizcd atomic frequency loa ■ The rate of spontaneous emission into the waveguide modes is given by 



71 = 2ngl A = ~^- 7o , (ID) 



where the second expression directly exhibits the scaling with the transverse extension of the waveguide. It is related 
to the atomic radiative cross section a a — 3X 2 /2n. A natural lower bound on the transverse mode size is at about 
A ~ (A/2) 2 (c.f. lowest order mode in box), implying a maximum achievable coupling ratio 71/70 ~ 1. In the range 
oa ~ A, one has a strong waveguide-atom coupling, which is manifested by the fact that the atom dissipates its 
energy equally into the waveguide and the free-space "lossy" modes. This situation could turn out to be a suitable 
basis to construct single-atom detectors or efficient light emitting diodes. 
For the following, it is convenient to introduce the "free pulse" operators 

/■oo 

9 uo a (t) = / dcug^aMe-'^-^ , (11) 







and the same for 60 (t), where ujq can be any frequency and will later be identified with the central frequency of the 
wave packet. The evolution of the atomic operators is given by the equations 



d 
~dt 



o- = -w- + 2a z9al0 (a (t - z A /c)e^ ZA ' c + b (t + z A /c)e-^' c ) e -*C«">-"A)t + |_ (12a) 



^-a z = -2 7 (c7 z + 1/2) 
dt 



9u 



{a (t - z A /c)e iu ° ZA / c + b Q (t + z A /c)e-**»A/ c )e- i ( Ua - u *') t + fc.c) + L , (12b) 



where 7 = 70 + 71 is the total decay rate of the atomic dipole. This is the generalization of the Bloch equations for 
the atomic dipole operators to describe the effect of a time-dependent external excitation in the case of a quantum 
driving field. For certain quantum states of the field, a closed set of equations can be derived from the above operator 
equations to calculate the matrix elements of the atomic operators. Note that only the free field appears in the 
dynamics of the atomic dipole operators. This is a direct consequence of the Markoff approximation, which implies 
that any change of the electromagnetic field by the atom cannot act back onto the atom. In other words, a photon 
emitted by the atom into the waveguide leaves the interaction region immediately and cannot be reabsorbed. There is 
no feedback mechanism as provided by the mirror in cavities. Hence the strong coupling between atom and waveguide 
yields a different dynamics with respect to the case of ordinary cavity QED. 

On inserting the solution of Eq. (12|a) into (g), the field operator can be obtained in the form 

M*) = M*o) + < cat W + ^ ack W + <f d W + C ise (t) , (13) 
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where the five terms represent respectively the free field, the field forward and backward scattered by the atom, the 
field radiated by an initially excited atom, and finally the vacuum noise which is coupled into the waveguide via the 
atomic dipole. They read 



a^ at (t) = g^ Q e- l ^-^ )zAjc \ dt' e l ^- UA+ ^ ]t ' / dt"e- i(w-WA+<7)t " ' 2a z {t")a a {t" - z A /c) 



(14a) 



,back / 



(*) = 9ugu e~ il * +Uo)ZA/c / dt'e i{ - u - UA+i ^ t ' / di"e- 4(w -^ +l7) *"2a 2 (i")&o(i" + ^/c), 



. rad / 



(t) = g^ WZA,c / dt'e t{ ^ A+ ^ ]t 'a^(t ), 



(t) = g u e- iuZA/c dt'e l ^- UA+ ^ ]t ' dt" e 7 *" 



(14b) 
(14c) 
(14d) 

We have set the initial time to — 0. Nevertheless, we still will use to sometimes in order to explicitly mark the initial 
value of the operators evolving in the Heisenberg-picture. In the double time integrals the order can be exchanged to 
carry out one of them. Substituting the solution in the electric field decomposition (^) , the integral over the frequency 
ui can also be performed. In this latter step, the lower bound of the frequency integration is extended to — oo. For 
ct > z — za, i.e. a light pulse has enough time to travel to the atom, one finds the following result for the waveguide 
electric field in the half space z > za 



hu>o 

AlTEcA 



, 1 OA 



ao(r) 

t+za/c 



2 A 10 



rft' e i(wo-^+n)(r+^/c-t') 2o - z ( t ') ao (t' _ ZA /c) 
+za/c 





dt > e i("0-UA+ l ' ( ){T + ZA/c-t') 2(7 ^ t ^ bo ( t , + Za I^ 



(15) 



Tiuja 
inecA 

+z A / 



2itgu>A e 



-(iuJA+l)(T + ZA/c) 



<T-{t ) + 



dt'e^'i-it') 



where r = t — z/c, and only the modes a w are considered as we are interested in the field outgoing in the +z direction 
(z > za)- In the subsequent lines one gets the free, the forward and the backward scattered, and then the directly 
radiated and the noise field. This is the central result that we use in the following to calculate measurable quantities 
in various experimental scenarios. One can draw a general conclusion already at this stage that the scattered terms 
are proportional, as naively expected, to the ratio <ta/A. 



III. SCATTERING OF A SINGLE LIGHT PULSE 



In this section we study the scattering process of a single light wave packet off a ground state atom. We calculate 
how the photons carried by the wave packet are redistributed among the forward, backward and transverse directions. 
The relevant quantity to be calculated is the expectation value of the Poynting vector defined as 

S(z,t) = Ai($ | fl (-)#) + £(-)B(+)|$ ) , (16) 
nuo Ho 

given in units of photon current in the forward and backward direction after the scattering event. In the present 
scalar model the magnetic field is just proportional to the electric field. The initial state l^o) is a simple product 
state of the waveguide field state, the atomic ground state, and the environment: 

l*o> = |^)|0 6 >|<7>|O e > , (17) 

The waveguide field describes a wave packet propagating in the directions +z, i.e. only the modes a u are excited 
with state \ip a ), while counter propagating modes b u are in vacuum state 1 0&) . As for the initial state \ip a ) we will 
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use Fourier-transform limited Gaussian pulses in a coherent state or in a single-photon Fock state, respectively. The 
mean Poynting vector of a free pulse is then 



SoM) = iV a 4Le-^ 2 (*-*/c) 2 



'2tt 



(18) 



where N a is the mean photon number carried by the wave packet and Q is the pulse bandwidth. 

As no field is radiated from an atom in the ground state, there is no contribution from the fourth term of Eq. (|l5|) to 
the Poynting vector. Similarly the free space environment as well as the backward propagating modes are in vacuum 
state. Therefore, all the terms containing the noise operator or b u (to) on the right-most side will vanish acting 
on the initial state \^o)- Thus the mean Poynting vector in the forward direction for z > za is formed from the 
superposition of the original light pulse and the one forward scattered by the atom, 



S(r) 



Q 

/27T 



7o 



( OA 



J 

JO 



u(t)\u(t)) + 

t+za/c 

dt'Re {T{t, t'){u(T)\2a z (t')\u(t' - z A /c))} 

> fT + Z A /c 

dt'T*{T,t') x 
dt"T{T,t"){u(t> - z A /c)\±a z (t')a z (t")\u(t" ~ z A /c)) 



(19) 



where r = t — z/c and 



) = cxp I — (t +t' - z a /c)(t + za/c - t') + i(uj - <jj a + h)(t + z A /c-t') 



(20) 



The Poynting vector is expressed in a product form with the original pulse shape function separated in a first term. 
To this end, we defined an auxiliary state 



\u(t)) = (27rfi 2 )-^e^ 2t2 ao WI*o) 



(21) 



As we work in the Heisenberg picture, the time dependence of this auxiliary state does not describe a dynamical 
process but follows from its definition. 

The reflected field in the domain z < za results from the light back-scattered by the atom (and additional quantum 
noise). As for the atomic radiation the two directions in the waveguide are equivalent, one gets an expression formally 
identical with the third term of Eq. (pj|) : 



^(back)^/^ _ 



n 
— 

/2tt 



'+z A / 
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PT +Z A /C r-T 

\ dt'T*{T\t') / 

Jo Jo 



'+z A /c 



dt"^(T , ,t")(u(t l - z A /c)\^ z {t')°z{t")\u(t" - za/c)) , (22) 



where r' = t — (2za — z)/c, expressing that the field propagates now in the —z direction. 



A. Coherent-state light pulse 



In the following we present how to evaluate the Poynting vector (19) in the case of wave packets initially in a 
coherent state. A definition of multimode coherent state pulses can be found in Ref. fll7| |. They have the property 
that they are eigenstates of the annihilation operators. A coherent-state pulse of bandwidth f2, centered initially at 
z = 0, can be defined through the eigenvalue equation 



a^(io)IV'a) = VK 



( 2 



-(w-w ) 2 /n 5 



(23) 
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where N a is the initial mean photon number in the wave packet. Its photon statistics is Poissonian around the mean 
N a with variance N a . 

It follows that |^o) is an eigenstate of the free pulse operator ao(t) and we have: 

\u(t)) = ^1*0) • (24) 
To evaluate the Poynting vector, Eq. (|l^), we need the one-time and the two-time averages of the population inversion 



operator. An equation of motion for (a z (t)) can be directly obtained by taking the mean of the Eqs. (12). Using 
again that |^o) is an eigenstate of the free-pulse operators, the usual form of the optical Bloch equations is found. 
We present it in the appendix. However, when the evolution of {a z {t)a z {t')) as a function of t is calculated from jl^), 
the operator a z (t') stands in between the free pulse operators and the state |^o)- We prove in the same appendix 
that ao(t — za/c) and tx(i') commute for t > t' and therefore the action of the pulse operator on its eigenstate can be 
easily carried out. This lemma allows the use of a form of the quantum regression theorem to calculate the two-time 
average {o z (t)o z (t')) . 

As it is shown in the appendix, an effective single-photon Rabi frequency g c s = g i ^ g (27rO, 2 ) 1 ^ 4 can be identified in 
the Bloch equations. This coupling constant suggests that the wave packet has an effective "volume" of \/2TTcA/fl. 
For larger bandwidths fl we get a stronger field per photon, which is a central quantity in cavity QED. However, 
the interaction times between the atom and the transform limited pulses get shorter in the same time, limiting the 
possibility of coherent operations on the atomic states. 

The photon number N a multiplies all the terms in the Poynting vector ([l9]). Beyond this simple linear dependence, 
the photon number has an influence on the evolution of the atomic population. Through this term, the atomic 
saturation introduces a nonlinear behavior into the Poynting vector, which can be a noticeable effect in the strong 
coupling regime even for relatively small photon numbers N a . 



B. Fock-state light pulse 



Now lets turn to a second example and study a wave packet of precisely one photon. The single-photon Fock-state 
is defined by 

|1«> = J^duj (^) 4 ^ ( ^ o)2/!)2 4(fe)|0»). (25) 

The definition in (||l]) leads to 

\u(t)) = \0 a )\0 b )\g)\0 e ) , (26) 

which, similarly to the coherent-state case Eq. (|24|), is time-independent. We now need the diagonal matrix elements 
(0, g\a z (t)\0, g) and (0, g\a z (t)a z (t')\0, g) where stands for both modes a and b and for the environment. Equation 
( |l2| ) implies that (0, g\a z (t)\0, g) = —1/2 for arbitrary time t. In fact, this value is the initial condition itself that 
makes the right hand side of the generalized Bloch equations vanish. Using again the commutation of the operators 
a (t — za/c) and a(t'), one can show that (0, g\a z (t)a z (t')\0, g) = 1/4 independently of i and t' . We get such simple 
solutions because we consider, formally, the Bloch equations when the atom is driven by a "vacuum pulse" . 

Having the exact solutions for the one-time and two-time averages, the integrals in Eq. ( [l9| ) can be analytically 
evaluated. The integration invokes the complex error function and the result itself is not very instructive. We omit 
this rather long formula here and only show typical results graphically in the figures of the next section. 



C. Single-atom transmittance and reflectance 



We continue the discussion of light pulse scatteringon a single atom by means of numerical examples. For this we 
evaluate the Poynting vector expressions (^) and ( p2] ) to determine how much of a light pulse is transmitted and 
reflected back by the atom. We introduce the transmittance, i.e., the transmitted mean energy (photon number) 
divided by the mean energy of the impinging wave packet. For an atom being on resonance with the pulse carrier 
frequency lja = wo we expect the strongest scattering. In the forthcoming examples the waveguide cross section is 
set to a a = -A. 

In Figure gthe single-atom transmittance is presented as a function of the pulse bandwidth on a logarithmic scale. 
The overall shape of the transmittance curve is very similar for a pulse initially in a coherent state with N a = 1 (solid 
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FIG. 1: Resonant single-atom transmittance as a function of the pulse bandwidth for pulses initially in coherent state with 
mean photon number N a = 1 (solid line) and in single-photon Fock state (dashed line) . The parameters are uja = (±>o (resonant 
pulse) and a a = A. 

line) or if it is in a single-photon Fock state (dashed line). The difference is manifest for bandwidths close to the 
transition linewidth, Q w 70. Apparently the effective saturation is reduced if exactly one photon is present instead 
of a distribution with mean one. 

As the incoming mean photon number is now 1 in this example, the transmittance gives directly the mean trans- 
mitted photon number. The atom presents an "obstacle" provided the pulse length is on the order of the natural 
lifetime or longer, i.e., for f2 < 70. The single-atom effect on the transmittance is approximately constant up to pulse 
lengths of few times the atomic lifetime. In this range the transmitted power is less than the half of the incoming 
one. The lower bound in the monochromatic limit is given by 

T=l-2-^+(^] , (27) 
7 V 7 / 

which is 4/9 « 0.44 in the example. There is an intermediate regime 70 < O < IO70 where the transmittance changes 
rapidly and it is noticeably smaller for a Fock-state pulse than for a coherent state. Finally, short light pulses, 
& 70, pass through the atom almost without being appreciably affected, realizing a sort of short pulse filter. A 
simple explanation is that a large spectral part of a broad-band pulse is far away from the atomic resonance, which 
reduces the effective coupling strength. Figure [l] therefore suggests that in situations where temporal resolution is 
required, for example in a detection scheme, the best choice for the bandwidth is f2 ss 70 where relatively short 
light pulses still experience significant attenuation. Longer pulses provide less temporal resolution without much 
improvement of the signal, whereas shorter pulses experience a significantly smaller change. 

Let us note briefly that the part of the impinging pulse which is not transmitted is either reflected or, most 
probably, scattered in the transverse lossy modes. The reflectance has an upper bound which is reached again in the 
monochromatic limit. It is (71/7) 2 — 1/9 ~ 0.11 in the numerical example considered for Fig. 1, that is at least with 
probability 4/9 the photon is lost from the waveguide. 

Due to saturation, there is a nonlinear effect in the transmittance. The more photons the light pulse carries, the 
more the atom gets excited. Population in the upper state of the atom reduces effectively the dipole strength, as 
follows from the term (a z ) in the Eq. which amounts therefore to a decrease in the amount of scattering. This 
effect is exhibited in Figure ^ where the transmittance is plotted against the mean photon number N a of the initial 
coherent-state pulse for various pulse bandwidths. Parameters are the same as in Fig. |IJ The starting points of the 
curves at N a = 1 correspond to four different points on the curve of the previous figure EL For longer pulses, see 
the curve (i) or (ii) in the figure, the saturation effect is small. This is because the effective coupling strength g e g 
decreases with the pulse length. Or, alternatively, one can view the same effect as a consequence of the photons 
arriving more distributed in time. For short pulses, on the other hand, the coupling constant becomes large enough 
to reach saturation even with weak incident light pulses. For example, curve (iii) corresponding to Q = 70 manifests 
a drastic nonlinear behavior even for photon numbers around 1. 

Figure H shows transmitted and reflected photon numbers for f2 = 70/IO (same as in curve (ii) of figure^) depending 
on the mean photon number of the pulse. The square root of the mean photon number is shown with error bars, 
representing the inherent quantum noise associated with coherent states. In order to detect the single-atom scatterer 
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FIG. 3: Transmitted (a) and reflected (b) number of photons as a function of the initial mean photon number N a - Error bars 
indicate the square root of the mean signal corresponding to the noise in a coherent state. The bandwidth is set Q — O.I70. 



with high probability, the change in the photon number (reflected or transmitted) must be larger than the noise. The 
figure reveals the minimum necessary photon number of the probe pulse to achieve the required operation regime. 
Photon numbers at about 20 in the probe pulse are sufficient to get a high quantum efficiency for the detector. 



D. Phase shift and pulse deformation 



When the carrier frequency is detuned from atomic resonance, the pulse can undergo a phase shift (dispersive 
scattering). This phase shift can be the basis of a single-atom birefringence in waveguides. This happens if the atomic 
dipole is parallel to one of the polarization modes and orthogonal to the other one. Only the first polarization mode 
experiences a phase shift whose magnitude will be estimated, based on Eq. (|l5|), in this section. 

Strictly speaking the study of the phase properties necessitates the definition of a phase operator and the investiga- 
tion of its time evolution. Here we will adopt a simplified treatment which allows perhaps a more instructive insight 
in the phase properties. Initially a coherent state pulse, defined in (|23|), is taken with the phase 0, i.e., the amplitude 
of all the components is real. We assume that the state stays almost a coherent state all along the evolution. This is 
fulfilled when the pulse induces little atomic population excitation. Then, formally, the inversion operator a z can be 
replaced by —1/2 times the unity operator, and the state |\&o) is indeed an eigenstate of the electric field operator, 
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Eq. (|l5|), at any time. The phase can be described by the argument of the corresponding complex eigenvalue, which 
we can write as 



£7W(r)|* )=i 



AnecA 



r a (<r)(l + %))l*o) 



(28) 



where ao(r) = y/N a (2nfl 2 ) 1 / 4 exp (— Jl 2 r 2 /4) describes the initial Gaussian pulse shape. The radiated and the 
backscattered terms, the third and the fourth terms in Eq. ([l5]), vanish because of the choice of the initial con- 
dition. The noise term is negligible as it is proportional to the atomic excitation, which is supposed to be small. 
Thus, the term h(r) stems exclusively from forward scattering. It can be generated in a simple form by changing the 
integration variable in Eq. ( |T5| ) as t' — > r + za/c — t' . One gets 



h(r) 



1 OA 

2 A 



■7o 



-+z A /c 



i)t' P -(7+2n 2 r)t' -in 2 *' 2 



(29) 



On inspecting Eqs. (Eq) and (|29|), two effects of the forward scattering on the coherent-state amplitude can be noticed. 
First of all, the coherent-state amplitude is shifted with a relative amount of h(r). Second, this shift depends on 
t = t — z/c which yields a deformation with respect to the initial Gaussian pulse shape. The upper integration bound 
in Eq. rt29) can be extended to infinity if the distance za between the initial pulse and the atom is very large compared 
to the pulse extension c/fl. Hence, the only dependence on r derives from the term 7 + 2fl 2 r in the exponent. As r 
varies within ±1/51 in the pulse, the ratio of 7 and H decides how much the pulse shape is distorted. 

For very long pulses (quasi-monochromatic excitation limit) the exponential decay with 7 dominates the integrand. 
The other exponentials can be neglected and the integral renders the well-known Lorentzian form of the single- 
atom susceptibility \ = (<ja/2A)/ ((luq — uja) + i"f)- Evidently, in this limit the result does not depend on r. For 
general relation between Q and 7, the integral has to be evaluated numerically. Let us first study the change of the 
coherent-state amplitude as a function of the detuning, and then its dependence on t. 

For the first case, we fix r = and plot the numerical solution of h(r = 0) in Figure |J. We show the real and the 
imaginary parts of h(0) as a function of the detuning loq — lua for two bandwidth values, Q = "fa and Q = O.I70. By 
analogy with the susceptibility, the real and imaginary parts can be regarded as an absorption and a dispersion-like 
curve, respectively. For reference, the Lorentzian function associated with the quasi-monochromatic excitation limit 
(lincwidth 7, oscillator strength (ta/2A) is also plotted with dashed lines. The three curves merge asymptotically for 
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FIG. 4: Change of the coherent-state amplitude as a function of the pulse carrier frequency detuning from the atomic resonance. 
Both the "absorption" Ke{/i(0)} and the "dispersion" Qm{h(0)} curves are plotted with solid lines for £1 = 70 and J! = O.I70. 
This latter corresponds to the larger h(0) values. Dashed curves show the Lorentzian shape for reference. 



very large detunings. The reason is that the frequency components of the pulse are more or less uniformly detuned 
from the atomic frequency, much like in the monochromatic limit. Hence, the phase shift can be well approximated 
from the analytic Lorentzian solution. Otherwise, for moderate detunings, the bandwidth dependence of the change 
in the coherent-state amplitude is quite apparent in the figure. The narrower the bandwidth, the larger phase shift is 
obtained. 

When the change in the coherent-state amplitude h(r) is small with respect to 1 and approximately proportional 
to i, that is, for ujq — lua ^ 7o 5 h directly gives the rotation angle in phase space, i. e., the phase shift. The phase 
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shift as a function of t is presented in Figure |5|. The detuning is chosen uiq — lua — 4,jo where the absorption is 
well reduced and the pulse power is almost completely transmitted through the atom. One can see in the figure that 
rotation angles in the range of 0.11 radw 6° are induced by a single atom for a a — A. Note that for il = O.I70 the 
phase shift is almost uniform, while it changes considerably for Q = 70. In both cases the deformation is asymmetric, 
showing that the front and the back of the pulse interact with an atom in a different state. 



IV. INTERFERENCE OF LIGHT PULSES 



One possible application of such a setup is nonlinear photonics, i.e. doing optics with few photons beyond beam 
splitters and mirrors Jl!| . As a prototype of this dynamics let us study now the "collision" of two pulses propagating 
in opposite directions. Instead of considering the modes 6 w (<o) initially in vacuum state, we will properly define an 
initial state that describes a wave packet centered at z = 2za- Then, this backward propagating wave packet and 
the forward propagating one encounter in the position of the atom at z — za where the atomic dipole can mediate 
an interaction between them. The calculations are based on the general result presented in Eq. (|l5|), which holds for 
any initial condition. 

We calculate the Poynting vector in the range z > za after the pulse collision. Here the field is composed of the 
initial pulse in the modes a u and its forward scattered part, superimposed with the back scattered part of the pulse 
in the modes b^. One gets an expression quite similar to the result in Eq. (|l9|), 



S{r) 



^1 _ 1 2 2 

-e 2 



/2tt 



u{t)\u{t)) 

t+za/c 



g 7 o J dt'2$le{F(T,t')(u(T)\2a z (t')\v(t' - za/c))} 

i-t+Za/c 

/ dt"T(r,t")(v(t' - z A /c)\A<j z {t')<j z {t")W - za/c)) 

Jo 



(30) 



just we needed to introduce a second state \v(t)), which is defined by 



\v(t)) = (2^ 2 )-M^ 2 (a (t)+e 2 ^/ c bo(t + 2z A /c)) |* ) 



(31) 



Let us consider first the case when the two counter-propagating pulses are in a coherent state. The initial state is 
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composed of the state \ip a ) as defined in Eq. (|23|), while the back-propagating modes are in the state \ipb) where 

i 

M*o)hM = V^*" 4 e -(—o) 2 /^ e ^ ZA /c^ b) (32) 

The last exponential term ensures that the position of the pulse is z — 2z^ at t = t . The state \v(t)) is proportional 
again to the initial state, 

\v(t)) = (v^ + e^y/Nt) l*o) • (33) 

One can immediately recognize that the appearance of the state \v(t)) in the second and third terms of Eq. (|30| ) 
gives rise to intriguing interferometric phenomena. Compared to the simple scattering scenario given by Eq. (|1S|), 
the second term, which is responsible for the absorption, can be enhanced. To this end, the counter-propagating 
pulses must have the same phase, i.e. (p = 0. Then the forward scattered pulse and the backscattered pulse interfere 
constructively yielding a reduced transmitted signal. On the other hand, the two pulses can interfere destructively 
when they black out each other in the position of the scatterer, i.e. for N a — Nf, and ip = 7r. The atom being in dark, 
this situation looks as if the atom were absent and the pulses propagate freely. This simple example demonstrates 
that there is indeed a non-trivial interaction between the pulses mediated by a single atom. Since the transverse 
scattering is missing, the outcome is definitely different from the coherent sum of the outcomes obtained with two 
single pulses. In between the extreme cases the outgoing mean photon number in the +z direction, denoted by N + 
varies as a function of the phase cp. This is shown in Fig. || where a difference of almost 80% is obtained for a 
bandwidth £1 = O.370 and a beam cross section <ja/ A = 1. 
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FIG. 6: Mean photon number N+ outgoing in the +z direction as a function of the relative phase <p of two identical counter 
propagating coherent-state pulses. The mean photon numbers are N a — Nt = 1, the bandwidth is £1 = O.370, and (ta/A = 1. 



Assume that we have a mean photon number N a fixed, e.g. N a = 1. Then the photon number N + outgoing in 
the forward +z direction can be tuned by varying the photon number N[,. This effect is represented in the Figure [?]. 
The well-pronounced minimum results from the concurrence of the second term of Eq. (|30|), which scales as \fN~t, and 
reduces the outgoing photon number, and the third term, which is linear in Nf, and increases the photon flux due to 
back scattered light. The manifest deviation from the linear behavior is due to atomic saturation. 

Let us finally consider the collision of two Fock-state pulses on the atom. The initial state is | \J/o) = l 1 a)|lfc)l.<?)|0 e ), 
where |lft) is defined by 

1 

\h) = J™ du f 4 e"(— "^/"V^/^t (to)|0 6 ) . (34) 

Initially there is no quantum correlation between the two Fock-state wave packets. As the phase in Fock states is 
completely undefined, one expects that the interference is missing in this case. This is true and can formally traced 
back to the fact that for the auxiliary state of Eq. ( p0[ ) , 



W)) = |0„,l t ) + |l o ,0i) , 



(35) 
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FIG. 7: Mean photon number N+ outgoing in the +z direction as a function of the mean photon number Nf, of the pulse 
incoming from the backward — z direction. Parameters are as in figure fj. 



the contribution of ao(t' — z^/c) ("forward scattered photon") and the contribution of b (t' + z^/c) ("back reflected 
photon" ) do not add up algebraically. Note that for coherent-state pulses this was different, the second "absorption" 
term of Eq. ([30]) got a factor \/N a + e lv y/Nb and the Bloch equations had to be solved only for the mean of atomic 
operators. Here, for Fock states, one has to solve a specific realization of the generalized Bloch equations (see in the 
Appendix). However, these equations are inhomogeneous and the solution is far from being the sum of two terms, one 
originating in |0 a , 1&) and the other in |l a , 0b). The numerically obtained result for the mean photon number outgoing 
in the +z direction is very similar, however, to the one obtained with one single impinging pulse. No interferometric 
enhancement or reduction of the outgoing photon number can be observed for pulses initially in Fock states. 

V. CONCLUSION 

For a field transversally confined in a tiny waveguide a single atom is able to have a significant effect on a light wave 
packet traveling across the atom. As limiting cases of the time-dependent ID scattering, one finds a transmittance 
reduction below 50%, strong nonlinearity in the transmittance even for very low energy pulses, and interferometric 
coupling of two wave packets with a visibility of up to 80%. All these effects clearly demonstrate that miniaturized 
waveguides or fibers coupled to an atomic dipole can be used as efficient single-atom detectors, or photon-photon 
couplers. As one might be able to control a single atom state on the quantum level, this could pave the way to genuine 
quantum devices for optics, as e.g. a quantum switch for light ^o|. Similarly one could envisage to construct single 
photon Bell state analyzers, as they would be needed for improved quantum cryptography or quantum teleportation 
setups. 

In our model we aimed to studying the fundamental nature of the interaction of a single atom with a waveguide field. 
A regime that can be referred to as "strong coupling regime" occurs when the transverse extension of the waveguide 
modes is close to the single-atom radiative cross section. The coherent interaction between the field and the atom 
can dominate damping. In contrast to cavity QED, however, the waveguide-atom coupling yields a dissipation-like 
evolution of the atom since the waveguide modes form a broadband continuum, similarly to the reservoir composed of 
the free-space radiation modes. On one hand, a single atom within the waveguide field becomes a significant scatterer, 
as featured by the effects presented in the paper. On the other hand, it is questionable if a weak quantum field can 
perform coherent population transfer in an atom. This is because a large effective single-atom Rabi frequency (g c s) 
is accompanied by the presence of a large damping rate (71) in the Bloch equations. In addition to this, the Rabi 
frequency depends on the interaction time defined by the pulse length. For example, for adiabatic passage techniques 
long interaction time is required, which inevitably leads to the reduction of the coupling constant. Our model equations 
provide a suitable ground to further study the dynamics of atoms coupled to one-dimensional continua of modes. 
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APPENDIX A: DERIVATION OF VARIOUS MATRIX ELEMENTS OF THE ATOMIC OPERATORS 

We prove in this appendix that the free pulse operator ao(t — za/c) commutes with the population inversion operator 
<J z (t') for times t > t'. This lemma allows the derivation of a closed set of linear differential equations for the required 
matrix elements of the population inversion operator and for that of the two-time product a z (t)a z (t'). We briefly 
present these equations which are formally very similar to the Bloch equations and the ones obtained by the quantum 
regression theorem in the case of driving an atomic dipole transition with a classical field. 

By definition, 



C=[ao(ti-z A /c),a,(t2)] = I dwe^-^-^ M*o),<r,(*2)] 

= J duje -^-^)(ti-zA/c) ^[aM,a z {t 2 )\ - gu) e~^ ZA ' c J^dt'e^-^' [<r_(tV,(t 2 )]} , (Al) 

where in the second step we used Eq. (||). Atomic and field operators, taken at the same time, commute, hence the 
first term vanishes. In the second term one can change the order of the integrals 

C = - e -^oZA/c J ' dt / e <(« - WA )f [ a _(t'),o z {t 2 )\ [ dug u e-* u - u °X tl - t ') . (A2) 



jo J 

The coupling constant g u being a slowly varying function of lu around the given optical frequency ujq can be taken 
out of the integral. The remaining inner integral amounts to a Dirac-delta function and 

C = -e-^A/c j" dt i^ -u A) t> [ a _tf) i<Tz fa)} 2irg Uo 5(t' - h) = 
Jo 

= -2ir 9ujQ e-^/ c e 1 ^-^ 0(t 2 - h) , (A3) 

where is the Heavyside step function, which proves our conjecture. 

The time dependent effect of the light pulse on the atomic dipole operators can be determined from the Eq. (|lj 
Depending on the initial state of the pulses, various matrix elements have to be calculated and inserted in the 
expressions for the Poynting vector Eqs. ([l9|), ( p2| ) and (^). In all the present cases, the required variables obey a 
set of linear differential equations that can be written in the form 

s = Bs + b . (A4) 

When the initial pulses are in coherent states, the evolution of the quantum mean of the population inversion 
operator is needed. It can be obtained by integrating the above differential equation with 



s 



3m {(*>+)} 



where 



-27 -2(V^Va + VNbCosip)g(t) 2^/Nlsm.ipg{t) 

2(y/Na + y/Nb cos tp)g(t) -7 to - lu a | . (,\:,) 

— 2y/Nb shi(pg(t) ll>a — ojq —7 
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Formally, this equation is equivalent with the Bloch equations obtained for an atom driven by a classical time- 
dependent excitation. The corresponding single-photon Rabi frequency reads 

9(t) = g eS e-^-^ 2 , (A7) 

with 5cff = 5LJ0 (2^ 2 ) 1 /4. 

The two-time averages in Eq. (n9h can be obtained by using the quantum regression theorem. For a fixed time t' 



prior to time t, the variables 

/ (*0> \ 

s(t)= (fiepJJ.tf)) , (A8) 
\(%m {a> + (t)}a z (t')) J 

obey the linear differential equation given in the Eq. (|Aj) with b T = (— j(a z (t')), 0, 0). 

When the counter-propagating pulses are initially in single-photon Fock-states, the matrix element (0 a , lf,|<7 z (i')|l) 
where |1) = |0 a , lb) + |l a ,0b) is needed. For symmetry reasons, the required element is just the half of (l|<7 z (t')|l), 



which can be obtained from the solution of (A4) with 



/ \ / -7 \ /-27 -4 5 (t) \ 

s(t) = 3?e{(0|«T^|l)} , b= 2g(t) , B= -7 u Q -u A . (A9) 
\Sm{{0\a'_\l)} J \ / \ lo a ~u - 7 J 

Note that non-diagonal matrix elements are involved in this set of equations. One can easily see that closed sets 
could be derived with increasing number of equations for higher photon number states, i.e., the presented approach 
is suitable to describe many other initial states as well. Finally, one gets for the necessary two-time average 

(l|*,(i)5-,(f)|l) \ /-7<lK(i')|l) 
s(*)= I (0\a'_(t)& z (t')\l) , b= g(t) 



(l\a> (t)a z (t')\0) / \ g(t) 



-27 -2g(t) -M*) 
B=( i{uj Q -uj A )-l | . (AlO) 

i(uj A - w ) - 7 
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